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ABSTRACT 

We investigate whether the recently-observed population of high-velocity white dwarfs 
can be derived from a population of binaries residing initially within the thin disk of 
the Galaxy. In particular we consider binaries where the primary is sufficiently massive 
to explode as a type II supernova. A large fraction of such binaries are broken up 
when the primary then explodes as a supernova owing to the combined effects of the 
mass loss from the primary and the kick received by the neutron star on its formation. 
For binaries where the primary evolves to fill its Roche lobe, mass transfer from the 
primary leads to the onset of a common envelope phase during which the secondary 
and the core of the primary spiral together as the envelope is ejected. Such binaries are 
the progenitors of X-ray binaries if they are not broken up when the primary explodes. 
For those systems which are broken up, a large number of the secondaries receive kick 
velocities ~ 100 — 200 kms -1 and subsequently evolve into white dwarfs. We compute 
trajectories within the Galactic potential for this population of stars and relate the 
birthrate of these stars over the entire Galaxy to those seen locally with high velocities 
relative to the LSR. We show that for a reasonable set of assumptions concerning the 
Galactic supernova rate and the binary population, our model produces a local number 
density of high- velocity white dwarfs compatible with that inferred from observations. 
We therefore propose that a population of white dwarfs originating in the thin disk 
may make a significant contribution to the observed population of high-velocity white 
dwarfs. 
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1 INTRODUCTION 

Oppenheimer et al (2001) have reported a population of 
high- velocity wide dwarfs in the solar neighbourhood which 
they claim to be halo objects. From their observations they 
infer a local number density of such objects to be ~ 10~ 4 
pc -3 much larger than previous estimates (Gould, Flynn & 
Bahcall 1998). Reid et al (2001) however, claim these stars 
are simply the high-velocity tail of the so-called thick disk 
population. This result is controversial; for example Hansen 
(2001) notes that the ages of the white dwarfs appear to 
have a spread, whereas the thick disk stars have similar ages 
(around 12 Gyr old). 

We investigate whether the observed population can be 
derived from the binary population within the thin disk. Re- 
cently, Koopmans & Blandford (2001) suggested that young 
white dwarfs could be ejected from the thin disk to the halo, 
through orbital instabilities in triple-star systems. In this 
paper, we consider a simpler scenario, where the secondary 
of a tight binary is ejected when the primary explodes as 
a supernova, breaking up the binary through the combined 
effects of mass loss in the supernova explosion and the kick 



received by the neutron star. To produce the high-velocity 
objects seen by Oppenheimer et al, one must impart kicks in 
the range 100 - 200 km/s to stars in the thin disk. As the kick 
velocity the secondary receives is similar to the star's binary 
orbital velocity at the moment of the supernova explosion, 
we require the stars to be in tight binaries, with separa- 
tions much smaller than the size of a pre-supernova red gi- 
ant. This in turn implies that the binaries must have passed 
through a common envelope phase. Once ejected, such sec- 
ondaries evolve into white dwarfs in a relatively short time 
if they are sufficiently massive. 



Our evolutionary scenario is discussed in Section 2 of 
this paper and described in more quantitative detail in Ap- 
pendix 1. In Section 3, we calculate the distribution of kicks 
received by the secondaries. In Section 4, we compute trajec- 
tories of secondaries within the Galaxy to determine whether 
they can produce a population of high-velocity objects sim- 
ilar to those seen by Oppenheimer et al. The likely produc- 
tion rate of such objects is discussed in Section 5 and Section 
6 gives our conclusions. 
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Figure 1. An evolutionary scenario for producing high-velocity 
main-sequence stars, some of which will later evolve into white 
dwarfs. Beginning with two main-sequence stars (phase 1), the 
primary evolves into a red giant (phase 2) filling its Roche lobe, 
leading to the onset of a common envelope phase (3) where the 
red-giant envelope engulfs the main-sequence star secondary and 
the helium-star core. The post-common-envelope system consists 
of a tight binary containing the helium-star core and the main- 
sequence star (phase 4). The compact star loses its helium enve- 
lope either via a wind or mass transfer, leaving a carbon-oxygen 
core (phase 5), which then explodes as a supernova, with the 
main-sequence star secondary often being ejected at high speed 
(phase 6). 



2 EVOLUTIONARY SCENARIO 

The evolutionary scenario for producing high- velocity stars 
considered in this paper is shown in Figure 1. We consider 
the evolution of binary stars where convective mass transfer 
occurs when the primary evolves off the main sequence and 
fills its Roche lobe. A common envelope phase follows. The 



secondary and the core of the primary spiral together as the 
envelope is ejected. After the common envelope phase, fur- 
ther mass transfer occurs if the primary (now a helium star) 
evolves to fill its Roche lobe for a second time. This (conser- 
vative) mass transfer increases the mass of the secondary but 
may also lead to the merger of the helium star and main- 
sequence star. Alternatively, if the primary fails to fill its 
Roche lobe, mass-loss occurs via a wind. In either case, the 
carbon-oxygen core of the primary remains. Providing the 
primary was initially sufficiently massive, it then explodes 
as a supernova. The secondary and the newly-formed neu- 
tron star become unbound in a large fraction of such bina- 
ries through the combined effects of mass-loss during the 
supernova explosion and the kick received by the neutron 
star on its formation. The kick velocity received by the sec- 
ondary will be similar to its orbital velocity at the moment 
of the explosion. Hence in very wide binaries, where a com- 
mon envelope phase is avoided, such velocities will be rel- 
atively low; too low to produce the high-velocity stars we 
require. However we will see in Section 4 that kick velocities 
<; 100 kms -1 are realised in a subset of binaries which have 
passed through a common envelope phase. A large fraction of 
the high- velocity secondaries subsequently evolve into white 
dwarfs, especially if a large fraction of secondaries are rel- 
atively massive. High-velocity stars should also be visible 
before they evolve into white dwarfs. For example Hooger- 
werf, de Bruijne & de Zeeuw (2001) and Magee et al (2001) 
both observe a population of high-velocity O and B-type 
main-sequence stars. 

It is possible that the helium star in the wider systems 
does not have sufficient time to lose its entire He-rich en- 
velope via a wind before exploding as a supernova. There 
would be two consequences of this. Firstly, the binary would 
be closer when the primary explodes as a supernova hence 
the secondary would have a higher orbital velocity. Secondly, 
the primary ejects more mass when it explodes, hence a 
larger fraction of binaries would be broken up. 

There are a number of restrictions on the initial sep- 
arations of binaries for the scenario described above. The 
system must initially be wide enough that the primary has 
evolved off the main sequence, and be convective, before it 
fills its Roche lobe. Yet the binary must not be too wide, 
otherwise no mass transfer will occur. We require that the 
main-sequence-star secondary not fill its Roche lobe at the 
end of the common envelope phase. If it does, we assume 
that the two stars merge and remove such binaries from fur- 
ther consideration. 

The evolutionary scenario described above is discussed 
more quantitatively in Appendix 1. 



3 POPULATION PRODUCED 

In this section we consider the population of objects pro- 
duced via the evolutionary scenario described in the previous 
section. If we assume that all stars of mass 8.3 Mq < M* < 
40 Mq produce type II supernovae, and that the IMF is 
given by dN/dM oc M -2 ' 35 , the majority of supernovae 
originate from stars less massive than 20 Mq. For illustra- 
tive purposes, we consider binaries of primary mass Mi = 12 
Mq . In Figure 2 we plot the initial separations as a function 
of secondary mass, M2, showing the various constraints for 
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Figure 2. Plot of initial separation a;, for M\ = 12 M@, as 
a function of secondary mass M2(both in solar units) showing 
constraints for the formation of a binary as described in Section 
3 and Appendix 1. 



the evolution of binaries as described in Appendix 1. The 
lowest of the three solid lines represents the minimum sepa- 
ration required to avoid mass transfer whilst the primary is 
still on the main sequence. The top two solid lines are the 
limits represented by equation (Al). Systems above the up- 
permost of the three lines will remain detached throughout 
the evolution of the primary. Those below the uppermost 
line, but above the middle line, will undergo mass transfer 
after the primary has evolved off the main sequence and be- 
come convective (and thus enter a common envelope phase) 
while those below the middle line but above the bottom 
line will undergo mass transfer while the primary is still ra- 
diative. This will also lead to a common envelope phase if 
M2 ;$ Mi/3. In this case, the common envelope phase will 
result in the merger of the primary core and the secondary. 
The dotted line is the limit set by requiring that the main- 
sequence star secondary remains detached immediately af- 
ter the common envelope phase (ie equation [A5], where to 
produce Figure 2 we have set oceAce = 0.5). The dashed 
line is the minimum initial separation to avoid the helium- 
star primary filling its Roche lobe when it evolves after the 
common envelope phase (ie equation [A6]). In the case of 
systems of primary mass, Mi = 12 Mq, mass transfer from 
the Helium star leads to a merger with the main-sequence 
star secondary for M2 Ss 3 Mq. Hence such systems will 
not produce high-velocity white dwarfs. The black dots in 
Figure 2 each represent systems evolving as described in the 
previous section, where the secondary is in a tight binary 
around the primary which explodes as a supernova. This 
region is limited from above by the requirement that the 
primary fills its Roche lobe (upper limit in equation [Al]) 
and from below by the requirement that the main-sequence 
star secondary remains detached immediately after the com- 
mon envelope phase (ie equation [A5]). For systems where 



Figure 3. Plot of kick velocities received by secondaries (in km/s) 
contained in binaries which arc broken up when the primary ex- 
plodes as a supernova. The binary population is that shown in 
Figure 2 where the initial primary mass is 12 Mq in all cases. 
The kick velocities are shown as a function of secondary mass M2 
(in solar units). 



the post-CE binary is sufficiently wide that the helium star 
does not fill its Roche lobe, this star loses its envelope in a 
wind, with the binary separation increasing slightly as given 
by equation (A10). In addition to the constraints shown in 
Figure 2, there is a lower limit on the mass of the secon- 
daries (M2 > O.8M0) if the secondaries are to evolve into 
white dwarfs during the age of the Galaxy. The dotted re- 
gion in Figure 2 will shift slightly for other values of oiceAce. 
However the range of binary separations at the moment of 
the supernovae explosion, and thus the distribution of kicks 
received by the secondary are similar. 

A similar picture to that shown in Figure 2 emerges for 
primary masses Mi ^ 12 Mq. For more massive primaries, 
the helium-star primary never fills its Roche lobe in the sys- 
tems of interest (as higher-mass stars produce higher-mass 
helium stars which evolve to have smaller radii). 

We now consider the effects of the supernova explo- 
sion on the binaries. We allow for the effect of mass-loss 
from the primary (we assume the post-supernova mass of 
the primary is 1.4 Mq). We also include the kick received 
by the neutron star at formation, using a kick distribution 
from Hansen & Phinney (1997). Approximately two-thirds 
of the binaries identified in Figure 2 are broken up. We show 
the kick velocity received by the secondaries, as a func- 
tion of secondary mass, in Figure 3. Here we have taken, 
as example, a Salpeter IMF for the secondary masses, ie 
dN/dM 2 oc M 2 ~ 2 ' 35 . We note from this figure that the kick 
velocity is often above 100 km/s. The distribution of kick ve- 
locities is similar for other values of primary mass. We also 
note that in all cases M2 <; 0.9 Mq, even though here we 
have assumed a Salpeter-like distribution for the secondary 
masses. The median mass for the secondaries is ~ 3 Mq 
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Figure 4. The positions of 1000 secondaries which have received 
kicks drawn randomly from the distribution shown in Figure 3 
(Mi = 12 Mq and initial secondary masses are drawn from a 
Salpctcr distribution). The trajectories have been followed for 10 
Gyr and the positions noted every 50 Myr. The positions are 
shown in cylindrical coordinates, with units of kpc. 



hence we would expect the vast majority of ejected secon- 
daries to have evolved into white dwarfs by today. In Section 
4 we use this kick distribution to compute the resulting spa- 
tial distribution of secondaries throughout the Galaxy. 



4 TRAJECTORIES OF HIGH-VELOCITY 
WHITE DWARFS 

Using the model for the Galactic potential described in Ap- 
pendix 2, we now compute the subsequent trajectories of 
the secondaries of the binaries which are broken up when 
the primary explodes as a supernova. For illustration wc 
again consider binaries of initial primary mass Mi = 12 
Mq, with the secondaries having a Salpeter-like distribu- 
tion of masses. Hence we use the kick distribution given in 
Figure 3. We will assume that the initial (circular) orbital 
radius of the binary in the Galaxy is drawn uniformly from 
the range 1 kpc - 10 kpc. We then apply to the secondary 
a kick drawn randomly from our kick distribution (and di- 
rected in a random direction with respect to the original 
motion within the Galaxy) and integrate its motion within 
the Galactic potential for 10 Gyr, writing out its position 
and velocity every 50 Myr. We consider 1000 trajectories. A 
plot of all the positions (in cylindrical coordinates) is shown 
in Figure 4. We note that a large fraction of positions are 
relatively close to their birth site. In Figure 5 we plot the 
tangential and radial velocities of all stars within the limits 
8 kpc < R < 9 kpc and —0.5 kpc < z < 0.5 kpc. We see that 
the distribution of velocities in the U — V plane is encour- 
agingly similar to that reported by Oppenheimer et al. The 
distribution shown in Figure 5 will clearly depend on the 



Figure 5. Plot of velocities (in kms -1 ) in the U-V plane for a 
subset of positions in Figure 4 given by 8.0 kpc < R < 9.0 kpc 
and —0.5 kpc < z < 0.5 kpc. 



distribution of kicks received by the secondaries, which will 
depend on the distribution of secondary masses as shown in 
Figure 3. 



5 PRODUCTION RATE 

In the previous section we showed that a reasonable set of as- 
sumptions for the kick distribution and the initial positions 
of the binaries yields a distribution of stellar velocities in the 
solar neighbourhood not unlike that seen by Oppenheimer 
et al. However, we have not yet addressed the issue of forma- 
tion rates of such systems. The formation rate is dependent 
on a number of uncertain factors: the global supernova rate 
within the Galaxy, the distribution of primary masses, the 
binary fraction, the distribution of secondary masses for a 
given primary mass, the distribution of initial binary sepa- 
rations, and the initial spatial distribution of binaries within 
the Galaxy. 

Given all these uncertainties, it is impossible to give a 
precise figure for the local number density of white dwarfs 
produced as described in this paper. Rather we can obtain a 
reasonable idea of the order of magnitude of the likely pop- 
ulation given a reasonable set of assumptions. Experimen- 
tally we note that some 0.3% of systems plotted in Figure 4 
are within the limits 8.4 kpc < R < 8.6 kpc and —0.1 kpc 
< 2 < 0.1 kpc. The volume of this region is V ~ 2 x 10 9 pc 3 . 
If we assume a (slightly high) supernova rate of 10 _1 /yr, the 
total number of type II supernovae over the entire life of the 
Galaxy is Asn ~ 1.5 x 10 9 . The local number density of 
ejected secondaries is given by Nsw/V ~2x 10~ 3 /2 pc~ 3 , 
where fa is the fraction of type II supernovae occurring in 
the sorts of binaries discussed in Section 3. In other words 
the process described here will produce the population of 
objects observed by Oppenheimer et al if f% J; 0.05. It is 
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quite possible that another set of assumptions would lead to 
an even higher production rate than the one derived here. 



6 CONCLUSIONS 

We have described a process producing a population of high- 
velocity stars from binaries in the thin disk. If the primary 
fills its Roche lobe after leaving the main sequence, con- 
vective mass transfer leads to a common envelope phase. 
The separation between the core of the primary and the 
secondary decreases significantly as the envelope is ejected. 
A second phase of mass transfer may occur if the primary 
(now a helium star) evolves to fill its Roche lobe. This phase 
is conservative and therefore increases the mass of the sec- 
ondary. Alternatively if the helium star fails to fill its Roche 
lobe, mass loss will occur via a wind. In both cases, the 
carbon-oxygen core of the primary remains. 

For initial primary masses, Mis > 8.3 Mq, this carbon- 
oxygen core then explodes as a type II supernova. Most of 
the binaries become unbound through the combined effects 
of mass-loss in the supernova explosion and the kick received 
by the neutron star at formation. 

We have computed the distribution of kick velocities for 
these secondaries. Considering a population of binaries ini- 
tially on circular orbits within the disk of the Galaxy, we 
apply kick velocities drawn from the distribution in random 
directions and integrate the motion of the ejected secon- 
daries within the Galaxy. We show that for a reasonable set 
of assumptions, we are able to obtain a population of high- 
velocity stars similar to that seen by Oppenheimer et al. 
There are number of uncertainties associated with calculat- 
ing a production rate of these stars, but for a reasonable set 
of assumptions, we show that this population may make a 
significant contribution to the population observed by Op- 
penheimer et al. 
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APPENDIX 1: EVOLUTIONARY SCENARIO 

In order for the late (convective) Case B mass transfer 
from the primary to occur, leading to the onset of a common 
envelope phase, we require that the primary fill its Roche 
lobe, but not until it has a radius larger than the maximum 
allowed for early (radiative) Case B mass transfer. This re- 
quires 



iW,rB(Mi,i) fl max ,B(Ml,i) 



/i(») 



/i(») 



(Al) 



where q Y — Mi^/A^i and -R ma x,B(Mi,i) is the maxi- 
mum radius allowed for Case B mass transfer for a star 
of mass Mi,i. The radius of star 1 when it fills its 
Roche lobe is given by Rij = fi(qi)ai, where fi(q) = 
0.49<? 2/3 /(0.6q 2/3 + ln(l + g 1 / 3 )) (Eggleton 1983). 

The mass of the primary after the common envelope 
phase is given by 



Mi 



M Hc (Mi,i) = oMi's 



(A2) 



where a = 0.125 and b — 1.4 (van den Heuvel 1994). 

The inspiral during the common envelope phase may 
be computed by equating the change in total energy of the 
two stars to the binding energy of the envelope up to an 
efficiency qce, ie E cnv — acE^E. Here E env can be written 
in the following form (Webbink 1984) 



GMi,i(Mi,i-Mi, E) 



(A3) 



A C E/l(l?i)ffli 

Combining equation (A2) with the expression for AE, 
after some rearrangement we arrive at the following expres- 
sion for the inspiral during the common envelope phase 



F = acs _ 2Mi,i(Mi,i-Mi, CE ) 

1_> ai acE Ace/2 (<ji)Mi,cEM2,i 



+ 



Ml,CE 



r 1 



(A4) 



We require that the main-sequence-star secondary not 
fill its Roche lobe at the end of the common envelope phase. 
If it does, we assume that the two stars merge and remove 
such binaries from further consideration. 



ai > 



-Rms(M 2j ce) 

-Fi^CE/2(<7CE) 



(A5) 



where Rms(M2,ce) is the radius of the main-sequence-star 
secondary and qce = Mi.ce/A^.cb- 

The primary will now be a helium star. Such stars will 
expand as they evolve and may, in turn, fill their Roche lobe 
if 



oce < 



fimax,He(Afl ! CE) 



(A6) 



/i(?ce) 

Under these conditions, conservative (Case BB) mass trans- 
fer will follow, where the envelope of the helium star will 
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be transferred to the secondary leaving the carbon-oxygen 
core of the primary. The masses of the two stars Mi t co and 
M 2 ,psn after such a phase of mass transfer are given by 



Mi,co =Mco(Mi,He) 



Mi 



= M 24 + (Mi,He - M 1;C o) 



(A7) 
(A8) 



The ratio of the separation after this phase to the sep- 
aration after the common envelope phase is given by 



fcE^F 



a PS N _ M ltCE M 2 



'ice 



M 2 M 2 



(A9) 



If, on the other hand, the He star never fills its Roche 
lobe, the envelope will be lost via a wind. Leaving behind 
the carbon-oxygen core, while the secondary mass will be 
unchanged (ie M 2i psn = A^ce). The ratio of the separation 
after this phase to the separation after the common envelope 
phase is given by 

_ apsN _ Mi,hc + M 2 ,ce 



F c 



acE Mi,co + M 2j psn 



(A10) 



The functions i? ma x,rB(M), i? max , B (M), i? max ,He(M), 
and Mco(M), were tabulated from Habets (1986), Paczyn- 
ski (1971) and Bressan et al (1993) and linear interpolation 
was used to evaluate them for a particular mass, M. 

In order for the primary to explode as a Type II su- 
pernova, we require that the core mass exceeds the Chan- 
drasekhar limit, ie Mi ; co > Men- For the assumptions here 
(equations [A2] and [A7]), this is equivalent to Mi,i > 8.3 
Mq. The mass-loss (from the supernova explosion) and the 
kick received by the neutron star may unbind the binary. 

If the binary is broken up, the kick velocity of the sec- 
ondary will be similar to its orbital velocity at the moment 
of the supernova explosion, which is given by 



V 2 ,PSN 



Mi 



•Wl.CO + M2,PSN 



G(Mi,cq + M 2 ,psn) 
apsN 



(All) 



$ = $ s + $d + $;, 



(A16) 



Following Paczynski (1990), we use the following choice of 
parameters: 

a s =0, b s = 0.277 kpc, M s = 1.12 x 10 10 M Q , (A17) 

a d = 3.7 kpc, b d = 0.20 kpc, M d = 8.07 x 1O 1O M , (A18) 



r c = 6.0 kpc, M c = 5.0 x 1O 1U M , 



(A19) 



Because of the cylindrical symmetry of the potential, the 
integration of the trajectories can be simplified to consider 
the evolution of the z and 7? components only, as given below 



dR dz 

It = VR ' Tt = r 



d,Vf 



dt 



dv. 



\dRj z 7?3' dt \dz) R 



(A20) 



where the z component of the angular momentum, j z 
Rv<t,. 



APPENDIX 2: THE GALACTIC POTENTIAL 



The Galactic potential can be modelled as the sum 
of three potentials. The spheroid and disc components are 
given by 

*« (*. z ) = ~^— r ° Ms ... m ( A12 ) 



® d {R,z) = 



(R 2 + [a s + (z 2 + b 2 y/ 2 } 2 ) 1/2 

GMd 

(R 2 + [a d + (z 2 + b 2 d y/ 2 ] 2 ) 1/2 



(A13) 



where R 2 = x 2 + y 2 . The component from the Galactic 
halo can be derived assuming a halo density distribution, 
Ph, given by 

Pc 



Ph 



(A14) 



1 + (r/r c ) 2 

where r 2 — x 2 + y + z 2 . The above density distribution 
yields the potential 



$> h 



GM C 



iln f 1 '" 



H atan 



(A15) 



where M c — \-Kp c f%. The total Galactic potential is the sum 



